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Abstract. The AdS/CFT correspondence has developed over the last years into a very useful and 
powerful tool for studying strongly coupled field theories at finite temperature and density. Of 
00 particular interest is the regime of near equilibrium real time evolution that can be captured via 

CN linear response theory. The AdS/CFT correspondence allows the calculation of retarded two point 

functions of gauge invariant operators by studying fluctuations around asymptotically AdS black 
Q holes. A major role is played by the poles of these holographic response functions: the quasinormal 

frequencies. I will review the appUcations of these ideas to the hydrodynamics of the strongly 
^ coupled quark gluon plasma and the holographic realization of strongly coupled superfluids. 
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O INTRODUCTION 
in 

en String theory is most often cited as the best candidate for the big TOE, the all encompass- 

CN ing theory of everything. Although there are many physicists who are quite outspoken 

^ sceptics of that claim, I think nobody can refute the fact that string theory is a verita- 

^ ble theory of everything in a slightly different sense: string theory unites and uses all 

^ theoretical concepts that are normally specific to certain sub-fields of physics, such as 

• r-j elementary particle physics, condensed matter physics and gravitational physics. 

r> One of these unifications so characteristic of string theory is the application of the 

AdS/CFT correspondence [1] to the realm of heavy ion collisions and to condensed 
matter physics. This sounds like a very exotic mixture of topics to talk about at a gravity 
meeting as was the ERE 2011. But the tools that are necessary to undertake a journey 
into the realm of heavy ion and condensed matter physics are indeed very well known 
to specialists working in general relativity. I will concentrate in my exposition on one 
particular of these tools: quasinormal modes of asymptotically anti-de Sitter black holes. 

Quasinormal modes arise because of the most basic properties of black holes: they 
swallow up everything leaving no traces back except for conserved charges. If we perturb 
a static asymptotically flat black hole space-time by any excitation, a scalar field, gauge 
field or metric fluctuation, these fluctuations can only end up being either swallowed 
by the black hole or radiated off to infinite. Quasinormal modes are the spectrum of 
linearized fluctuations with outgoing boundary conditions at spatial infinity and infalling 
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FIGURE 1. Quasinormal modes are small fluctuations obeying infalling boundary conditions on the 
horizon, (a) In asymptotically flat space times they also obey outgoing boundary conditions at infinity, (b) 
In asymptotically AdS space times the boundary conditions at r — > oo can be either Dirichlet, Neuman or 
mixed ones. In the context of Holography the they are fixed to be the poles of the holographic retarded 
correlators 



boundary conditions on the black hole horizon. The mathematical problem one has to 
solve is one dearest to the hart of every theoretical physicist: an eigenvalue problem 
for a differential operator. However, contrary to the eigenvalue problems one encounters 
normally, the particular boundary conditions make the the differential operator a non- 
hermitian one and therefore the eigen frequencies turn out to be complex numbers of the 
form 

COn = n„ - iTn . (1) 

The real part is the frequency with which the perturbation oscillates whereas the imag- 
inary part is the rate at which this oscillation is damped. The fluctuations are damped 
if r„ > 0, otherwise an instability in form of an exponential growth in time appears. 
In asymptotically AdS spaces spacelike infinity acts like a reflecting boundary and Neu- 
mann, Dirichlet or some form of mixed (Robin) boundary conditions have to by supplied 
there. Quasinormal modes for asymptotically AdS black holes have first been consid- 
ered in [2] and in the context of the AdS/CFT correspondence in [3, 4]. The holographic 
dictionary singles out a specific boundary condition corresponding to the poles of the 
holographic Green's functions [5, 6]. 

Unfortunately the mathematical theory of quasinormal modes lacks far behind its her- 
mitian cousin. Few general statements about the spectrum (1) can therefore be made, e.g. 
only for some simple equations stability can be proved. The quasinormal mode spectrum 
in most cases has to be calculated using some numerical approximation method. For the 
AdS/CFT applications it is however enough to study the properties of some of the lowest 
lying quasinormal modes. It turns out that these are intimately connected to the hydrody- 
namic behavior of strongly coupled field theories. Very useful reviews on quasinormal 
modes are [7, 8] 



THE ADS/CFT CORRESPONDENCE 



Before discussing the role of quasinormal modes we need first to understand the basics 
of the AdS/CFT correspondence. We will be very brief in summarizing the main ideas. 
The fundamental statement is "Type IIB string theory on AdSs x S5 is dual to ,JV — 4 
supersymmetric gauge theory". Let's try to understand this statement in a crash-course 
style. The jy = 4 supersymmetric gauge theory is a non-abelian, four dimensional 
quantum field theory whose field content consists of six scalars, four Majorana fermions 
and a vector field. They all transform under the adjoint representation of the gauge 
group SU (N). It features four supersymmetries and this fixes all the couplings between 
the different fields. As it is a gauge theory physical observables are gauge invariant 
operators such as tr(F^v^'^^)- The global symmetry group S0{6) acts on the scalars and 
the fermions (in the SU{4) spin representation of 5(9(6)). Because of the high amount 
of supersymmetries the theory also has exact conformal symmetry on the quantum 
level and this leads to four additional conformal supersymmetries , such that we have 
2 X 4 X 4 = 32 supercharges in total. 

The dual theory is a theory of gravity (this is what type IIB string theory is) but living 
in quite a few more dimensions, ten as opposed to the four the field theory knows about. 
But five of these ten are easily got rid off: the isometries of the ^5 part of the geometry 
form S0{6). The S5 is the geometric realization of what appears as an internal, global 
symmetry group in the field theory. Type IIB string theory has two supersymmetries in 
ten dimensions, they are Majorana- Weyl spinors of same chirality with 16 components 
each, which gives the 32 supercharges necessary to match the field theory counting. 

The field theory is characterized by two parameters, the gauge coupling gYM and 
the rank of the gauge group A'^. The dual string theory has a string coupling g.^ (the 
amphtude for a string to split in two) and a fundamental length scale Is, the string scale. 
Furthermore the geometry is determined by a scale L determining the curvature of the 
AdSs as R= — 20/L^. The AdS/CFT correspondence relates these parameters in the 
following way: 

sImN - ^ (2) 

- gs (3) 

The AdS/CFT correspondence is therefore a strong weak coupling duality: for weak 
curvature we have large L and therefore also large 't-Hooft coupling. In this regime 
of weak curvature stringy effects are negligible and we can approximate the string 
theory by type IIB supergravity. If we furthermore take the rank of the gauge group 
A'^ to be very large we can also neglect quantum loop effects and end up with classical 
supergravity! This is the form of the correspondence most useful for the applications to 
many body physics: classical (super)gravity on{d+l) dimensional Anti-de Sitter space 
is the infinite coupling and infinite rank limit of a gauge theory in d dimensions. Here we 
have allowed ourselves to be already a bit more general. Once we have understood the 
original example based on the maximally supersymmetric four dimensional field theory 
we can conjecture that every theory of gravity plus some suitably chosen matter fields 
on (d+l) dimensional Anti-de Sitter space is a dual to a certain quantum field theory 



in d dimensions. In fact we might even be a bit more brave and delete the words "dual 
to a" in the previous phrase. This is the point of view taken in the applications of the 
AdS/CFT correspondence to the world of solid state physics. The additional matter fields 
are then chosen to reflect a particular symmetry content of the underlying quantum field 
theoretical system one is interested in. Having this in mind we will forget from now on 
some of the non-essential ingredients, such as supersymmetry and the extra dimensions 
in form of the S^. 

For the applications to quantum field theory the most useful foliation of Anti-de Sitter 
space is the so called Poincare patch in which the line element takes the form 

ds^ = -zr^i-dt^ + dx^) + ^dr^. (4) 

The space on which the dual quantum field theory lives is recovered by taking the limit 
dsqpj = limr^ac r^^ds^. This is why sometimes it is said that the dual field theory lives 
on the "boundary" of AdS and why the AdS/CFT correspondence is also referred to as 
"Holography". The radial coordinate has however a physical interpretation as an energy 
scale. We can identify the high-energy UV limit in the field theory with the r — > oo limit 
in the AdS geometry, whereas the low-energy IR limit is r — )> 0. 
The asymptotic behavior of the fields in AdS has the form 

$ = r-^- (^)o(x) + 0(r-2)) + r-^+ {^^ (x) + 0{r-^)) . (5) 

The exponents A± obey A_ < A_|_ and depend on the nature of the field, e.g. for a scalar 
field of mass m they are A± = ^{d ± \/~(P^+~AnPl?) . 

We now would have to evaluate the path integral over the fields in AdS keeping the 
boundary values 4>o(jc) fixed. The result is a functional depending on the boundary fields 
4>o(x). Now the boundary field ^o(-^) is interpreted as the source J{x) that couples to a 
(gauge invariant) operator &{x) of conformal dimension A+ in the field theory 

Z[y]= / exp(-j5[4)]) (6) 

Connected Green's functions of gauge invariant operators in the quantum field theory 
can now be generated by functional differentiation with respect to the sources 

In the limit in which the gravity theory becomes classical, i.e. the large and large 
coupling Sym^ P^th integral is dominated by the classical solutions to the field 

equations and logZ can be replaced by the classical action evaluated on a solution of the 
field equations. In this case the coefficient 4>i (x) of the asymptotic expansion (5) is the 
vacuum expectation value of the operator sourced by 4>o. 

(^W)oc4>i(;c). (8) 

To explicitly compute ^i{x) we need to supply a second boundary condition, so 
far we have fixed only the asymptotic value <I>o. For time independent solutions we 



demand regularity in the interior of the (possibly only asymptotically) AdS space. For 
time dependent solutions we need a different boundary condition. The metric (4) has a 
(degenerate) horizon at r = and it was argued in [9] that for time dependent solutions 
the appropriate boundary conditions are infalling ones. In particular it was argued 
that the Green-functions that are calculated with these infalling boundary conditions 
correspond to retarded ones in the field theory, e.g. on the level of two point functions 
the retarded correlator is given by 

Gr{t,x) = -i®{t) ([^(?,f),^(0,0)]) . (9) 

That does not yet lead to quasinormal modes as long as we are in pure AdS space. 
But that changes as soon as we consider a black hole with planar horizon topology, i.e. 
a black brane in AdS. Its line element is 

ds^/L" = r2(-/(r)dr2+^je2) + ^— ^, (10) 

m = 1-^, (11) 

where for simplicity of notation we have rescaled t and x. This metric as a non degenerate 
horizon atr — rn — TtT with planar topology. The Hawking temperature is T. A natural 
conjecture is that this metric is now the dual to the strongly coupled field theory in 
a state of thermal equilibrium at temperature T. More precisely it corresponds to the 
plasma phase of the dual gauge theory [10, 11]. For a conformal theory in Minkowski 
space that any temperature puts the theory in the plasma phase. Because of conformal 
symmetry every finite temperature is as good as any other. This is reflected in the gravity 
background by the fact that we can rescale the coordinates r,t,x and to set T = 1. 

Suppose now that ^» is a small fluctuation around the black hole background. We 
expand it in plane boundary waves 

^ir,t,x) = I ^Mco:k)e-'^'+^'^-'F^^kir) . (12) 

For every fixed (O^k the linearized equation of motion for the fluctuation boils down 
then to an ordinary second order differential equation for F^-^. In turns out that in 
asymptotically (non-extremal) AdS black hole backgrounds this equation is of Fuchsian 
type, i.e. having at most a finite number of regular singular points in the complexified 
r-plane. The point at infinity is a regular singular point with characteristic exponents 
A±. Imposing infalling boundary conditions ~ g-i(o{t+r^) j-j^g horizon (we use 
a tortoise coordinate here dr^ — drjf{r) such that the horizon sits at — oo) the 
asymptotic expansion of a solution takes the form 

- A((0,X) r-^- [1 + (9(l/r)] +fi((0,X) r-^+ [1 + 0(l/r)] . (13) 

The retarded two point correlator turns out then to be given by the ratio of the two 
expansion coefficients A, 5 [9] 

G,(a,:k)=K?^. (14) 



This leads to the following holographic definition of quasinormal modes in AdS: they 
are solutions of the linearized fluctuation equations that correspond to poles of the 
holographic retarded Green's function in the complexified frequency plane. 

Retarded two point functions are the central objects in linear response theory. If we 
perturb a system by switching on a small perturbation represented by a source J{(oJc) 
that couples to the operator the response is given by 

If G/; is a holographic Green's function computed from an asymptotically AdS black 
hole background with regular horizon it has indeed only poles in the complex frequency 
plane and the integral can be evaluated according to Cauchy's theorem 

m,x)) = @{t)Y^R,,{k)j{co„{k))e-^''"'+^'^^ (16) 

n 

where Rn are the residues of Gr at the poles a)„. As long as all quasinormal frequencies 
lie in the lower half of the complex w-plane the response is decaying exponentially fast. 
If one (or more) modes happen to lie in the upper half plane we find a mode whose 
amplitude grows exponentially with time, getting soon out of the linearized regime and 
indicating an instability. 
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FIGURE 2. The ten lowest lying quasinormal modes for a scalar field excitation (from [12]). 

If the operator ^ is a conserved charge then the spectrum of quasinormal modes 
contains special modes whose dispersion relation fulfills 

limQ}H(k)=0. (17) 

The reflect that fact that the total charge can not dissipate away. We call these special 
modes hydrodynamic. These modes play a distinguished role in the applications of the 
AdS/CFT correspondence. 




FIGURE 3. Nature confronts us sometimes with complicated objects such as the one in the upper left 
corner. In simple situations, as for example in collisions at ultrarelativistic energies, it is sufficient to 
replace it by the simpler object in the right upper corner. The same philosophy can be applied to the 
strongly coupled quark gluon plasma! 

THE HOLOGRAPHIC QUARK GLUON PLASMA 

Let me come now to the first application of the ideas outlined in the previous section. 
What can quasinormal modes tell us about the physics of the quark gluon plasma? 

Strong nuclear interactions are on a fundamental level described by QCD, a non- 
abelian gauge theory with gauge group SU{3) and three quarks transforming in the 
fundamental representation. It is therefore a theory of the same type as the ,yV = 4 
supersymmetric gauge theory, however QCD is of course not conformal. Its gauge cou- 
pling runs leading to asymptotic freedom in the UV and confinement in the IR. Its 
spectrum consists only composite color neutral objects such a Mesons and Baryons 
at low energies. If we are able however to increase the energy density to values such 
that \/£ ~ Aqcd we expect the theory to undergo a phase transition (or at least a fast 
crossover) to a deconfined regime. Such high energy densities can be generated in exper- 
iments colliding heavy ions, i.e. gold or lead atoms completely strip off their electrons, at 
relativistic energies. Two such experiments are currently running, a dedicated one at the 
Relativistic Heavy Ion Collider (RHIC) in the USA and in Europe at the Large Hadron 
Collider (LHC) which every year dedicates its november run to heavy ion collision. 

One of the spectacular findings of these experiments is that the droplet of hot matter 
created in the collision behaves as an almost ideal fluid. It viscosity to entropy ratio 
seems to be the smallest one ever observed for a fluid. The value is nowadays most often 
cited in units of J=, e.g. in [13] we find 



l<4;r- <2.5. 

s 



(18) 



These values are obtained by assuming that the quark gluon plasma can be described by 
hydrodynamics. Simulations based on relativistic hydrodynamics where rf/s is an input 
parameter can be used to compute the particle multiplicities as function of transverse 
momentum. Fitting these to the data allows the extraction of a value for rj/s. This is 
of course not the only input parameter, initial conditions play a role as do the time 
from when on the hydrodynamic description is supposed to be valid. The latter is often 
referred to a thermalization time scale. Data can be well fitted assuming thermalization 
times 0.6/m/c < T < Ifm/c [13]. 

The temperatures reached in these heavy ion collisions are about two to three times 
the deconfinement temperature. This is a rather low value and indicates that QCD is still 
a rather strongly coupled system at these scales. Perturbative estimates of rj/s and the 
thermalization time are in fact to high to be compatible with experimental data. We seem 
to be in need of a genuine strong coupling description of the quark gluon plasma. Since 
numerical methods based on lattice QCD are mostly confined to Euclidean signature 
they are of limited help in calculating transport coefficients such as the viscosity. 

We do have however an accessible model for a strongly coupled non-abelian gauge 
theory in the plasma phase via the AdS/CFT correspondence. It is of course quite 
different from QCD in its details but at least it might serve as a toy model. Let us see 
if we can subject this toy model to some simple tests that probe its applicability to the 
quark gluon plasma created in heavy ion collisions. 



Hydrodynamics 

Let us briefly recall the basics of hydrodynamics. It can be understood of an effective 
field theory. The equations of motion are simply the conservation laws of energy- 
momentum and possible some conserved currents. The energy-momentum tensor is 
expanded in the number of derivatives that act on the thermodynamic variables such 
as pressure p, energy density e and the fluid four velocity u^. It is described by the 
equations 

= d^T^'"', (19) 
T^^ = (e + p)MM„v^^^Mv + ^Mv^ (20) 

T^^ = -T](7^^-CP^%mP, (21) 

(22) 

with the transverse projector P^^ — u^u^ -\- T]^^ and the shear tensor a^"^ ~ 
P^^P"^^ {daUp + dpUa — ^Vap^p^'^)- Assuming the validity of the equation of state 
p = p{e) at each space time point, i.e. local thermal equilibrium, we have four dynami- 
cal equations for the four variables p and the normalized u^. At this level there are two 
transport coefficients describing dissipation, the shear viscosity rj and the bulk viscosity 
^. It is not too difficult to see that they have two eigenmodes with dispersion relations 

©=-/— (23) 

e+p 

(0 = ±vsk-rsk^. (24) 



They are the shear mode and the sound mode, = §f is the speed of sound and Fg = 
2(e'+p)' sound attenuation constant. Using the thermodynamic relation e + p = sT 
we see that 77/5 determines the diffusion constant of the shear mode. 



Quasinormal modes 

Let me concentrate on the shear mode (23). It obviously fulfills (17). If the AdS black 
brane is a strongly coupled plasma then we should be able to find such a mode in its 
quasinormal mode spectrum. This is indeed the case. Since the source of the energy 
momentum tensor is the metric we switch on a metric perturbation gMN ^ gMN + ^MAf- 
We can chose the three momentum to lie in the z-direction. The relevant perturbation 
can then be packaged into a gauge invariant combination Z = ^^^i + which in the 

coordinate u — ^^-^ fulfills the equation [14, 15] 

Z"{u) + fiV/il \ 2f^2 " + ?f^^ " = , (25) 

uf{u) [q^fW - w^) uf{uY 

where {w,q) = 27tT{co,k) are dimensionless. In these coordinates quasinormal modes 
fulfill infalling boundary conditions at the horizon at w = 1 and Dirichlet boundary 
conditions z = at m = 0. 

Shear viscosity to entropy ratio: The low lying spectrum is depicted in figure (4). 
It shows the imaginary parts of the four lowest lying quasinormal frequencies in the 
shear channel. There is one hydrodynamic one, whose frequency goes to zero with the 
momentum whereas the others are gapped. The hydrodynamic mode is purely imaginary. 
The gapped modes also have non-vanishing real parts (not plotted). Hydrodynamics 
applies in the low frequency long wavelength limit. In this limit 0), ^ — > the equation 
(25) can be solved analytically and one finds [16, 17] 

co=-i—k^. (26) 

The shear viscosity to entropy ratio of the strongly coupled ^ = 4 plasma in the large 
limit is therefore 

s An 

This we can take as a first success of our toy model! As mentioned this value is 
indeed compatible with the experimental estimation and in any case much better than 
perturbative estimates. In fact not only the ^ = 4 theory gives this value but many 
holographic models with less symmetries, e.g. less supersymmetry or non-conformal 
backgrounds give precisely the same value. There is a certain amount of universality 
to this result. So far all holographic models with isotropic horizon and at most two 
derivative actions give this value [18]. For some time it was conjectured that l/(4;r) is 
a fundamental lower bound [19]. If one includes however higher derivative terms the 
value can be even smaller [20, 21] (see [22] for a review of the effects higher derivative 





0.2 




J h 1 








(} 




41.1 




4 

:u> 




A 
2Ii 




2 


1 


l.h 




1 




FIGURE 4. Quasinormal modes in the shear channel as a function of momentum. In the lower part 
of the figure the imaginary parts of the lowest lying modes are plotted. The ungapped mode shows a 
diffusion like dispersion relation at low momentum. For at higher momentum or shorter wavelengths it 
crosses with the first gapped mode. The at these short wavelengths the system does not behave according 
the hydrodynamics anymore. The upper part shows the residue of the diffusive mode. It goes to zero as 
the momentum is increased. 



terns have on r\/s). Very recently it has been shown that certain anisotropic backgrounds 
can violate the conjectured bound even at the two derivative level [23]. 

Validity of hydrodynamic approximation: However from the figure we see that the 
system does not behave in a hydrodynamic way down to arbitrary small length scales. 
This poses another test on our toy model. The length scale at which this happens should 
be much smaller than the size of the fireball created in heavy ion collisions. Otherwise 
the toy model would predict that the fireball can not be described by hydrodynamics! 
The fireball can be estimated to have a diameter of 10/m. The temperature lies between 
200 MeV and 300 MeV The value at which the hydrodynamic mode crosses the first 
non-hydrodynamic mode is ^ = 1.3(2;rr). We use the conversion factor 1 fm^^ = 197 
MeV and find a length scale of / = (1.3 * 2 * :;r * 200/197)^1 ^ 10"^ fm. This is good 
news! The fireball is approximately 100 times bigger than the length scale at which our 
strongly coupled toy model shows non-hydrodynamic behavior. Our toy model predicts 
therefore that hydrodynamic simulations should be good approximations to capture the 
flow of the quark gluon plasma. 

Thermalization time scale: In full generality this is of course a very complicated 
problem and relies on a good modelling of the regime far from equilibrium beyond 



the linear response theory accessible to the quasinormal mode analysis. Whatever this 
complicated far from equilibrium dynamics might be we can be sure however that at 
sufficiently late times the system equilibrates in form of quasinormal mode excitations. 
Hydrodynamics is the late stage of this process in which the lowest, hydrodynamic 
mode dominates the time development. We can estimate this time by comparing the 
first gapped with the hydrodynamic mode and demand that the response is dominated 
by the hydrodynamic mode from the time^ Xh on. From (16) we can estimate 

\0g\RHJ{0}H)\-\0g\RlJ{0ix)\ 

Tff = :^ (^o) 

i W -i 1 

This depends of course on the particular form of the source J. By a judicial choice 
of source this time can be made as long as one wishes. For a variety of reasonable 
sources investigated in [24, 15] it turned out however that Xh ^ 0.14 fm/c. This again is 
good news. If this time had turned out to be much larger than 1 fm/c without having to 
adjust the sources in a special way it had meant that our toy model predicted generically 
thermalization times larger than the values inferred from experiments. 

I also want to comment on one more feature of the quasinormal mode spectrum 
depicted in the figure. If we go to large momenta the prolongation of the hydrodynamic 
mode grows faster and faster. In fact it seems that the ratio 0)/k grows as well and goes 
way beyond the causality bound of 1. This looks a bit worrisome at first. But if we 
look at the residue of this mode we see that it goes to zero. So in the high momentiun 
regime this mode effectively ceases to exist. Furthermore to study the issue of causality 
one should really study the front velocity in which one keeps the frequency real and 
searches for poles in the complex k plane [25]. The large frequency behavior of the 
mode hydrodynamic mode is perfectly causal with a front velocity Vf = 1 (see figure 
(5). 

Let us summarize what we have seen from this simple considerations. The strongly 
coupled c/T = 4 gauge theory seems to be a reasonable good toy model for the strongly 
coupled quark gluon plasma. It predicts that hydrodynamics is applicable to a quark 
gluon plasmaball as small as 10 fm and the quasinormal modes show decay rates that 
at least allow for a fast thermalization. It also predicts a shear viscosity to entropy ratio 
compatible with the values inferred from data via hydrodynamic simulations. These are 
some of the most basic features of the holographic model of the quark gluon plasma. 
There exist of course many more instances where quasinormal modes play a crucial role 
in holographic models of the quark gluon plasma. A comprehensive review on AdS/CFT 
applications to the physics of the strongly coupled quark gluon plasma can found in [26]. 



HOLOGRAPHIC SUPERFLUIDS 

Superfluidity is one of the most emblematic phenomenon of many body physics. It is 
therefore quite natural to ask if a holographic implementation of superfluidity can be 
found. On a fundamental level it is the question of spontaneous symmetry breaking and 



This is sometimes also called the scrambhng time. 
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FIGURE 5. Keeping the frequency real one can solve for complex momentum values. At small frequen- 
cies the behavior is the typical diffusive one 0) — —iDk^ whereas for large frequencies the front velocity 
limftj^oo j^lp^ = 1 can be read off. 



it has been answered to the affirmative in [27, 28]. A simple model of a holographic 
superfluid is the one presented in [29]. The model was first formulated in a 2 + 1 di- 
mensional setting, i.e. in a four dimensional asymptotically AdS space. This is however 
not crucial and 3 + 1 dimensional models can easily be constructed. The quasinormal 
mode spectrum of this model has been investigated in [30] and in the following we will 
briefly outline the construction of the model and the main features of its quasinormal 
mode spectrum. 

One characteristic of superfluids is the phenomenon of second sound. Superfluid He- 
lium II for example has the capacity of transporting heat not in the usual diffusive form 
described by a heat equation but through propagating waves. This is the phenomenon of 
second sound. We will review how this can be investigated in a holographic superfluid 
with the help of quasinormal modes. 

We start with a four dimensional AdS black brane with metric 

ds^ = -f{r)dt^ + T7T + [dx^ + dy^) ■ (29) 
/('") 

where f{r) = j2 — This geometry has a regular horizon at vh = M^/^L^/^ with a 

Hawking temperature of T = 

In addition we will consider a simple Abelian Higgs model on this background with 
Lagrangian 

^ = - IfmnF^^ - - Dm^D^^' , (30) 

Fmn is the field strength of an abelian gauge field Am and ^ is a complex scalar field 
charged under the t/(l) gauge symmetry. Its mass is chosen to be = —2L^ which 
would be tachyonic in flat space but is perfectly stable in asymptotically AdS spaces. 
We will also consider only the so called decoupling limit in which we neglect the back 



reaction of the gauge and scalar field onto the geometry. This can be justified in the limit 
of large charge of the scalar field [29] . 

Now one can look for non-trivial solutions to the field equations in which the temporal 
component of the gauge field Aq = 4> and the scalar field take the following form 
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(31) 
(32) 



According to the holographic dictionary ji is interpreted as a chemical potential, n the 
charge density, as the source for a dimension 2 operator i//2 = k^^ijiux {^i)- 

The underlaying conformal symmetry allows to fix one of the dimensionful parame- 
ters, we can chose e.g. /i = 1 and view the variation of the dimensionless combination 
/i = 1^ ^ as variations in the temperature. 

It turns out that for high temperatures there is a unique solution with trivial scalar field 
^ = and gauge field 4> = /i (l - ^) . 

Below a certain critical temperature Tc — 0.587/i there is however a second solution 
with non-trivial scalar field profile with y/i = but ^f2 ^ 0. This is the sign of sponta- 
neous symmetry breaking: although there is no source for the operator ^2 it does have a 
non trivial vacuum expectation value. 
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FIGURE 6. The figure on the left shows how the lowest quasinormal modes in the scalar field channel 
move as the temperature is varied. At the critical temperature the lowest two modes sit precisely at the 
origin. There they contribute to the hydrodynamics. For even lower temperatures the move into the upper 
half plane indicating an instability towards condensation the the scalar operator The right figure shows 
the condensate as a function of Temperature. In this broken phase there are no unstable modes. 

In the normal phase the the quasinormal mode spectrum decomposes into three 
non-interacting channels: the scalar field channel, the transverse gauge field and the 
longitudinal gauge field channel. The transverse gauge field channel shows a gapped 
spectrum, the longitudinal one features one ungapped hydrodynamic mode. It reflects 
charge diffusion and has a diffusive behavior of to = —iDk^ with diffusion constant 
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. The most interesting behavior can be found in the scalar field channel. For 



high temperatures the spectrum is gapped but the gap diminishes as the temperature is 
lowered. At the critical temperature two modes becomes ungapped, one for the scalar 
field ^ and one from the complex conjugate field ^. If one lowers the temperature even 
further but stays in the normal phase with vanishing expectation value for the charged 
scalar operator these modes move into the upper half of the complex frequency plane. 
They are therefore exponentially growing modes that signal an instability towards the 
formation of a scalar condensate! 

In the broken phase with non trivial ^ background the scalar field channel and the 
longitudinal vector field channel mix. Dividing the scalar in its real a and imaginary r\ 
part the equations of motion become 

(33) 

. II f ri 2/\ , 2 0)^ 1^^> 2i(0(p 

= fat" + —a/-('^+m>Aat-^a^-2i(0'¥7]-4^'(j>G, (35) 
P VP / P 

. + -2*^).,+ f., + 2,W,. (36, 

subject to the constraint ju/ + ^«x' = 2i{\f/' rj — YV') ■ construct three solu- 

tions obeying infalling boundary conditions at the horizon and fulfilling the constraint. 
A fourth one can be obtained by taking a pure gauge field configuration. These span a 
basis of the space of solutions to (33) which can be assembled in a 4 x 4 matrix. The 
quasinormal mode conditions is now that the determinant of this matrix evaluated at the 
boundary of the asymptotically AdS space vanishes. 

The goldstone theorem suggests that now there should be a massless mode in the 
spectrum with dispersion relation 

co = ±Vsk-irsk^ , (37) 

V.V is the speed of second sound and Fs is its attenuation constant. Careful search of the 
determinant condition shows indeed that this mode is in the spectrum and that the purely 
imaginary mode that represented diffusion in the unbroken phase is now gapped 

co^ -iy-iDk^ . (38) 

So the only hydrodynamic modes in the broken phase are the Goldstone modes corre- 
sponding to second sound. The square of the speed of sound and the gap in the purely 
imaginary mode close to the critical temperature behave as 



The speed of sound and the attenuation constant can be obtained numerically by fitting 
the dispersion relation (37) to the poles of the determinant as shown in figure (7). As can 
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FIGURE 7. The square of the speed of sound and its attenuation constant as determined from fitting the 
lowed quasinormal mode to the dispersion relation (37). 








FIGURE 8. The higher, non-hydrodynamic modes move continuously but not smooth through the 
second order phase transition. 



be seen the attenuation constant rises steeply close to the critical temperature but stays 
finite as one approaches it. 

The counting of hydrodynamic modes is as follows: in the normal phase there is one 
diffusive mode, precisely at the critical temperature there are three modes, the diffusive 
one and the two scalarmodes that become ungapped at the onset of the phase transition. 
In the broken phase the purely imaginary mode is gapped and there are two ungapped 
modes corresponding to second sound. 

Finally it is also of interest to see how the non-hydrodynamic modes behave as one 
goes through the phase transition, this is depicted in figure (8). The second and third 
scalar modes are shown as the critical temperature is approached. They behave in a 
continuos but non-smooth way. For temperatures lower than the critical one the modes 
do not anymore approach the real axis buy move roughly parallel to it as one lowers the 
temperature further. This suggests that no new instabilities arise. 

Applications of the AdS/CFT correspondence to condensed matter physics are re- 
viewed in [31, 32, 33, 34, 35] 
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